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Abstract -Charge carriers in single and multilayered graphene systems behave as chiral particles 
due to the particular lattice symmetry of the crystal. We show that the interplay between the 
meta-material properties of graphene multilayers and the pseudospinorial properties of the charge 
carriers result in the occurrence of Klein tunneling and cloaking for rhombohedral stacked multi- 
layers. We derive an algebraic formula predicting the angles at which these phenomena occur and 
support this with numerical calculations for systems up to four layers. We present a decomposi- 
tion of an arbitrarily stacked multilayer into pseudospin doublets that have the same properties 
as rhombohedral systems with a lower number of layers. 



Introduction. — The electronic properties of the one 
atom thick graphene crystal has been the subject of sev- 
eral recent papers jl] since its experimental isolation. [2] 
Not only is the crystal a promising candidate for semi- 
conductor physics, the electron behavior also mimics that 
of Dirac particles and can therefore be seen as an inter- 
esting table top realization of a two dimensional quantum 
relativistic system. The Klein paradox, a unit transmis- 
sion probability through potential barriers of any height 
or width, was one of the first characterizing phenomena 
from QED predicted and subsequently observed ex- 
perimentally i4j. But also optical properties as Fabry- 
Perot resonances [5] and the negative refraction index 
that makes graphene a metamaterial [6] are remarkable 
properties of the crystal. The stacking of two layers of 
graphene, the so called bilayer graphene (BLG), although 
being only weakly bound, changes fundamentally the elec- 
tronic properties. The Klein paradox as observed in mono- 
layer graphene (MLG) is replaced by the cloaking [7] of 
hole states inside a potential barrier resulting in a suppres- 
sion of transmission which is sometimes called anti-Klein 
tunneling. For trilayer graphene (TLG), Klein tunneling is 
present if the layers are orthorhombic stacked. For Bernal 
stacking however Klein tunneling is absent j8l|9]. 

The occurrence of these tunneling phenomena is the 
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consequence of the lattice induced chiral nature of the 
charge carriers. When the Fermi energy of the electrons is 
low enough, in SLG the dispersion of the electrons in the 
vicinity of one of the Dirac points K is linear in recipro- 
cal space. This allows for the introduction of pseudospin, 
which is the lattice induced analogue of conventional spin 
from the Dirac theory. Extending this concept to BLG, 
the dispersion near the Dirac point can be approximated 
as being parabolic but due to symmetry arguments still 
leads to a spinorial Hamiltonian describing the electrons 
as chiral particles with a pseudospin [TUIITT] . 

Some recent papers have discussed the electronic struc- 
ture |121I13| and known concepts such as trigonal warping 
and the Berry phase [TIHTB] of graphene multilayers. In 
this paper we generalize the discussions of Katsnelson et 
al. and Gu et al. [7] to an arbitrary number of layers 
and to an arbitrary stacking order for a PN junction. We 
find that the low energy behavior can be expressed as a 
system of non interacting pseudospin doublets, each with 
a specific chirality and derive a simple algebraic expression 
for the angles at which Klein tunneling and cloaking can 
be expected. 

Structure of the multilayer. — A system of n lay- 
ers of graphene can be stacked using a multitude of differ- 
ent stacking sequences. Graphene consists of two trigonal 
sublattices, called a and /3, and it therefore suffices to con- 
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sider only the relative position of these sublattices of the 
different the layers. Due to the periodicity of the crystal, 
there are only three ways a layer can be placed with re- 
spect to the bottom one which wc label as A, B and C. 
It is possible to place it directly above the bottom layer 
(A), to shift it once (B) or twice (C) by the interatomic 
distance in the direction of the vector between the a and 
fi atoms. Bernal stacking (ABA) [TTlfTS] is found to be the 
most stable combination. Rhombohedral stacking (ABC) 
is however also possible and is observed for a small number 
of layers [I9ll20j . In between these two stacking structures, 
each different combination leads to a different electronic 
structure. To find the low energy electronic spectrum of 
the different possible stackings for an arbitrary number of 
layers, one can use the decomposition method as described 
by Min et al. [H]. Notice however that they do not dif- 
ferentiate between structures that are symmetric under in 
plane mirroring and thus represent the same system. The 
total number TV (n) of physically different stacking possi- 
bilities for n > 3 layers is given by 
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where [aj refers to the nearest integer smaller than a. 
The second term of this expression takes into account the 
mirror symmetry of a stacking sequence and is absent in 
the discussion of Ref. [ [H]]. In FigUwe show for n = 4 
the 3 stacking possibilities together with their respective 
energy spectrum near the K point. 

Rhombohedral multilayers. For n rhombohedral 
(ABC) stacked graphene layers, the spectrum consists of 
two bands that touch at the K point and 2n — 2 bands that 
are located at higher energies. If we only include the near- 
est neighbor interlayer hopping, so only hopping between 
the j3i and a^+i sublattices, the effective Hamiltonian of 
such system near the K point is given by a 2n x 2n matrix 
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with a ~ [ux, cFy) a vector of Pauli matrices, vp ~ lO^m/s 
the Fermi velocity in single layer graphene, k the wave 
vector and r is given by 



1 

hvp 




71 



(3) 



where [21j 71 = 377meV is the interlayer hopping param- 
eter. For this kind of stacking it is possible to introduce 
a two band low energy approximation which yields the 



Fig. 1: (Colour online) (left) Electronic energy spectrum of the 
three physically different stacking possibilities of the graphene 
tetralayer with only nearest neighbor interlayer hopping. The 
energy is expressed in units of the interlayer hopping energy 71 
and the wave vector is expressed in units of a~^, the inverse of 
the nearest neighbor interatomic distance, (right) Schematic 
representation of the relative positions of the sublattices a and 
/3 of each layer for the stacking orders shown on the right. The 
yellow line corresponds to the interlayer hopping taken into 
account in the calculations. 



Hamiltonian 
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, (4) 
(5) 



where (j)k = arctan (ky/kx) is the angle of the wave vector 
k with the normal chosen perpendicular to the PN junc- 
tion and ax,y are the components of the pseudospin associ- 
ated with this two dimensional Hamiltonian which are the 
respective Pauli matrices. The validity of this approxima- 
tion is shown in Fig. [5] for n = 2, 3, 4, 5. In this figure we 
show the dispersion relation obtained by the Hamiltonian 
in Eq. ^ which consists of 2n bands. Superimposed in 
this figure we have plotted the dispersion relation from the 
Hamiltonian in Eq. ([5]) as dashed curves. The two band 
spectrum is in good agreement with the 2n band spectrum 
for low energy and near the K point. 

The two band Hamiltonian in Eq. ([5|) describes a chiral 
particle for which the pseudospin spins n times as fast as 
the wave vector k. Notice that for the angles 
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the sine or the cosine in Eq. 



,0,...,n-l}, (6) 
([5|) vanish for m even or 
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Fig. 2: (Colour online) Energy spectrum near the Dirac point 
of rhombohedral stacked graphene multilayers up to five layers. 
The black solid curves correspond to the tight binding energy 
with only nearest neighbor and next to nearest neighbor hop- 
ping and the red dashed curves correspond to the two band 
approximat ion . 



odd respectively. At these angles the Hamiltonian com- 
mutes respectively with Ox or Oy making them conserved 
quantities. 

When electrons impinge on a PN junction with an angle 
of incidence given by Eq. (I6|), conservation of pseudospin 
allows the electron to be reflected only if the reflected state 
has the same pseudospin as the incident state. The angle 
of the reflected state is given by 0^ = tf — i^^. an therefore 
the wave vector of the reflected electron must rotate by 
an angle of Ac/) = tf — 2(/)fc. Since the pseudospin rotates n 
times as fast as the wave vector, we must have riA0 = Z27r, 
/ G N, in order that reflection is allowed by pseudospin 
conservation. Using Eq. (|6]) yields that the pseudospin of 
the reflected state is parallel with that of the incident state 
for angles (^"^ when the difference n — m is even, while the 
pseudospin is opposite when n — m is odd. 

When the Fermi energy (E) of the incident electron is 
less than the potential step (y) of the PN junction, the 
sign of the wave vector of the propagating hole state inside 
the junction is opposite to the sign of the incident electron 
wave vector. This is the result of charge conservation at 
the steps' edge and gives rise to a negative refraction index 
as found for monolayer graphene [5l[6] making it a meta- 
material. Due to the change in sign of the wave vector, 
the angle of the wave inside the potential region also flips 
sign, making the pseudospin inside the junction to rotate 
in the opposite direction. This is schematically illustrated 
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Fig. 3: (Colour online) Schematic representation of the pseu- 
dospin matching at a PN junction for (a) bilayer and (b) tri- 
layer graphene with Fermi energy E = F/2. The red arrows 
show the rotation of the pseudospin with respect to the in- 
cident angle in the reciprocal plane outside (left) and inside 
(right) the junction. The blue arrows indicate the angles for 
which the conservation of pseudospin results in Klein tunneling 
(arrows directed to the right) or cloaking (arrows directed to 
the left). 



in Fig. [3] In this case the angle of refraction is given by 



arctan 



kx [Th) 



(7) 



where (T^) is the wave vector of the hole state and 
Th = \E — V\ corresponds to the kinetic energy of the 
hole. Due to electron-hole symmetry, when E ~ V/2, the 
refractive angle is exactly opposite to the angle of inci- 
dence. Following a similar argument as before, one finds 
that for angles (/)™ the pseudospin is opposite to that of 
the incident state when n — ?7i is even, but it is the same 
when n — m is odd. 

A mismatch in the pseudospin inside and outside the po- 
tential step was invoked earlier to explain Klein tunneling 
and cloaking in monolayer and bilayer graphene [3]l7] . Fol- 
lowing the above argument, one can conclude that for an n 
layered rhombohedral stacked system, cloaking is present 
at angles given by Eq. (j6]) when n — m is even, while Klein 
tunneling is present for angles given by Eq. ^ when n — m 
is odd. Table [T] lists the special angles for n up to 5. 

For arbitrary stacking, Min et al. [12] showed that the 
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Table 1: Angles for Klein tunneling (KT) and cloaking (CI) for 
multilayer graphene with n = 1, . . . , 5 layers. 
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Hamiltonian can be decomposed in a set of independent 
pseudospin doublets of the form 

Hji ^ fc/' [cos {Ji(j)k, ) cTx + sin ( Ji(/)fe, ) ay] , (8) 

where ki is the wave vector corresponding to the propagat- 
ing low energy band of the i"* pseudospin doublet and Ji 
is the number of layers taking part in the doublet which 
corresponds to the chirality of this doublet. The num- 
ber of pseudospin doublets Nd depends on the details of 
the stacking, but the sum of the chiralities equals the to- 
tal number of layers in the system. In this way one can 
decompose the low energy structure of any multilayercd 
system in a set of Njj non interacting doublets. There- 
fore, there arc Njj modes of propagation for low energy. 
As long as the symmetry of the system remains intact, it 
is not possible to scatter between different modes. The 
Hamiltonian given in Eq. is that of a system of Ji 
rhombohedral stacked layers of graphene and because the 
different modes don't interact, the occurrence of cloaking 
and Klein tunneling is the same as before for a multilayer 
with n = Ji. 



Transmission probability. — The two band Hamil- 
tonian in Eq. ([5]) has a plane wave solution given by two 
propagating waves, one right and one left moving, and 
2n — 2 evanescent waves. The wave vectors of these plane 
waves are the solutions of the equation 



(9) 



where e = E / (hvp)" and ky is the transverse 

wave vector. The solution of this equation is zLkj with n 
different values for kj. The plane wave solution can be 
written as a two spinor 
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where the latter is a matrix formulation of the spinor with 
matrices 
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and 



1 • • • 1 ; J 
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To find the transmission probability for a PN junction, one 
has to equate the plane wave solutions and all the deriva- 
tives up to n — 1*'' order of the region before the junction 
(region I) with those of the region behind it (region II) at 
the junction's edge at a; = 0. This leads to a set of n two 
component equations: 



ViEiCi = PhShCii 



(13) 
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where the matrix £ is evaluated at a; = 0. Normalizing 
the incident wave on the right propagating wave before 
the junction by putting j ~ 1 and applying boundary 
conditions ajj = and a^jj ~ for j ^ 1 to suppress the 
non normalizable plane wave functions, the transmission 
(T) and the reflection probability (i?) are given by 
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The numerical results for the transmission probablity for 
multilayers with n = 1 up to 4 are depicted in Fig. U] 
as function of the energy of the incident electron and the 
incident angle. The expected angles for Klein tunneling 
and cloaking are confirmed by our calculations. 

Conclusions and remarks. — The combination of 
the chiral nature of the charge carriers and their meta- 
material properties induce Klein tunneling and cloaking at 
specific angles for rhombohedral stacked graphene multi- 
layers. For an arbitrary stacking sequence, the low energy 
behavior of the electrons can be decomposed in indepen- 
dent chiral doublets with a chirality J that act as if it's a 
rhombohedral multilayer with n = J layers. In this way, 
for any arbitrary stacking sequence, one can predict the 
occurrence of Klein tunneling and cloaking. Note how- 
ever that we limited ourselves to nearest neighbour hop- 
pings and neglected other less important hoppings that are 
present in a real multilayer. The latter results in e.g. like 
trigonal warping |16j that will effect our results for very 
small energies. Furthermore, the use of the two band ap- 
proximation limits the energy range to about SOOmeV. At 
high energies, additional modes of propagation need to be 
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Fig. 4: (Colour online) Transmission probability through a PN junction in multilayer graphene for n = 1, . . . , 4 layers. The top 
row shows the angle dependend transmission for electrons with E = V/2. The bottom row gives countourplots of the angle and 
energy dependence of the transmission. 



taken into account, changing the transmission properties 
for high junctions and high Fermi energy 
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